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ABSTRACT A method for predicting the noise field of a propfan operating 
at angle of attack to the oncoming flow is presented. The method takes 
advantage of the high-blade-count of the advanced propeller designs to 
provide an accurate and efficient formula for predicting their noise field. The 
formula, which is written in terms of the Airy function and its derivative, 
provides a very attractive alternative to the use of numerical integration. 
A preliminary comparison shows rather favorable agreement between the 
predictions from the present method and the experimental data. 


1.1 Introduction 

The need to limit both in-flight cabin noise and community noise has 
prompted the development of methods for accurate prediction of acous- 
tic performance of propfans. A particular aim of these methods has been 
the prediction of the noise field of a propfan operating at angle of attack to 
the oncoming flow. Experimental evidence [1] suggest that the noise radia- 
tion characteristics of a propfan at angle of attack are significantly different 
from those of a propfan at zero angle of attack. These differences stem from 
the changes which occur both in the aerodynamics and aeroacoustics of a 
propeller whose axis is inclined with respect to the direction of oncoming 
flow. 

One of the most important of these changes is the azimuthal variation 
of aerodynamic loading experienced by a propeller at angle of attack. Mea- 
surements [2], as well as numerical simulations [3], have shown that the 
loading variation or “unsteadiness” can be quite substantial over one period 
of revolution. Physically, the modulation occurs because propeller blades 
experience periodic variations in the local flow incidence. Acoustically, the 
loading unsteadiness creates a noise field which is axially non- symmetric. 
Methods which predict this effect have been in use for some time (see Han- 
son [4], for example). 

Another feature of a propeller operating at angle of attack, which is ab- 
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sent at zero angle of attack case, is the cross-flow convective phase effect. 
This effect, first pointed out by Mani [5], may be explained as follows. For 
an inclined propeller, the oncoming flow has a nonzero velocity component 
in the plane of rotation. As a result, an acoustic source point fixed to the 
surface of the propeller “sees” a medium whose convective velocity varies 
periodically. This means that the convective phase of the sound radiated 
from this point is also periodically modulated. Krejsa [6], recently devel- 
oped a prediction scheme which accounts for this effect by assuming that 
propeller sources are acoustically compact in the chordwise direction. 

In this paper, a method which incorporates both the loading unsteadi- 
ness and cross- flow phase effects in a straightforward manner, without the 
simplifying assumptions such as chordwise compactness, is presented. This 
method, which is based on a moving-medium variant of the Ffowcs Williams 
and Hawkings [7] equation, relies on a frequency-domain formulation simi- 
lar to that employed by Hawkings and Lawson [8] or Hanson [9], however, 
it differs from them in that the resulting Fourier integrals are evaluated 
asymptotically in the ‘large-blade-count limit” rather than in the far- or 
near-field limits. The advantage of this approach lies in the fact that it 
provides a uniform representation for the sound field in both the near- and 
far-fields while, at the same time, it affords the flexibility of using realis- 
tic blade geometries and aerodynamic loading distributions without their 
associated computational penalties. 


1.2 Analysis 

The starting point for the analysis is the moving-medium solution of the 
Ffowcs Williams and Hawkings equation for the acoustic pressure, p(x, <), 
as given by Goldstein (TOj; 


/ ^ f + °° f . , w [ + °° f . dG , . . , 

p(x,t) = - / pov n ——ds(y)dr+ / /n<^— i«(y) rfr 

J-oo Js(t) Dt j - CO Js(t) vyi 


Do 

Dr 




d 

d Vi ’ 


(la) 

(16) 


where p 0 is the ambient density, v n the normal component of the surface 
relative velocity, / the amplitude of the aerodynamic loading, and Tij the 
Lighthill stress tensor. S(r) and V(r) represent the propeller blade surfaces 
and the volume surrounding the blades, respectively. njVare components 
of the outward unit normal to 5. It should be noted that the convention for 
the direction of the unit normal is opposite that of Goldstein’s. The medium 
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convection velocity is denoted by Uo. y and x are the source and observer 
coordinates, respectively, r denotes the source time and t the observer (i.e., 
retarded) time. The three terms in Eq. (la) represent contributions from 
the thickness (monopole), loading (dipole) and stress (quadrupole) sources 
of noise. In this paper, however, the contribution from the last term will be 
ignored. The extension of the analysis to include quadrupole sources will 
be discussed elsewhere. 

G = G(x,f/y,r) in Eq. (la) denotes the free-space moving-medium 
Green’s function which may be written in the following form; 

G=^~6(t- T - gc R/ c o), (2a) 

0e(r) = i (* ~ MoJ , «(r) = (M 0 2 b + Pi) 1 ' 2 , 

P 0 

M 0 r (t) = M 0i e it e >( T ) = - Pa = 1 _ M oi> 

Mot = R(t) = |x - y(r)|. (26) 

Here Mo R is the component of the medium convection Mach number in the 
radiation direction, e^’s the components of the unit vector in the radiation 
direction, R the length of the radiation vector, and Co the medium speed 
of sound. Parameters g c and « represent the effects of medium convection 
on the retarded time and the spherical spreading rate, respectively. The 
explicit dependence of various parameters on the source time, r, is indicated 
where necessary. 

In order to simplify the description of the source motion, the 1-axis is 
chosen to coincide with the propeller shaft (see Fig. (1)). As a result the 
motion of the sources is confined to the transverse planes described by the 
2— and 3— axes. In this coordinate system, the medium convection Mach 
number can be written as: 


Mqi — M 0 cos a, M 02 = 0, M 03 = M 0 sin a, (3) 

where a is the propeller angle of attack with respect to the oncoming flow 
as shown in Fig. (1). 

Given the blade geometry and aerodynamic loading, the acoustic pres- 
sure, p(x,f), may be computed from Eq. (la) for any observer location. 
This, of course, is the so called time-domain approach (see Farassat [11], 
for example) which entails solving transcendental equations for the source 
time, r. Alternatively, one could use a frequency-domain approach which 
involves expanding the pressure field, p(x, t ), in terms of its complex Fourier 
harmonic components, pjt(x), i.e., 
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(4a) 


+ OO 

p(x,t)= Y, Pk(*)*~ iktU > 

k= — oo 

with the individual harmonic components given by: 
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(46) 


where fl is the propeller angular speed. In writing Eq. (4b), the derivatives 
of G with respect to the source coordinates from Eq. (la) have been re- 
placed by its derivatives with respect to the observer coordinates through 
the use of the simple relations: 


dG _ dG D 0 G _ DpG 

dyi dxi } Dr Dt 

Using the usual phase relationship arguments it can be shown that for 
B identical blades only pC s for which k — mB (where m is an integer) 
contribute to the infinite sum in Eq. (4a). This contribution is simply B 
times the contribution of a single blade. Therefore, from here on B will 
explicitly appear in the expressions for p m Q and, correspondingly, S will 
denote the surface of a single blade. 

To further simplify Eq. (4b), the spatial derivatives, d/dx{ , of G may be 
rewritten in terms of the temporal derivative, d/dt , through the use of the 
chain rule; 


dG _ _ [_L d(g c R) d_ l _ dR 1 dtC 

dxi Cp dxi dt + R dxi k 3x,- 

The temporal derivatives can now be removed using integration by parts 
and the integral over t } which involves a delta function, can then be easily 
evaluated to yield, 



r2*/0 r 

b(x) = / / Qmonopo\e(x,y,T) e ,mBn( - T+ScR /c»') ds(y) dr 

J 0 Js(t) 


/*2*/0 r 

/ / £?dipole(x, y, r) e imBCl(r +gc R/C a ) dfi(y) 

Jo JS(t) 


( 7 ) 


where Q m0 no P oie and Qdipoie represent the expressions for the “amplitudes” 
of the two acoustic source. Note that the limits on the integration over r 
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have been changed to reflect the fact that an interval of size 2w/Q in t is 
mapped exactly into an interval of size 2 ir/Cl in r. The expressions for the 
monopole and dipole source amplitudes are given by: 


Mb, +$«)], (8a) 

9c^0i) + 

-^(Mo R M 0 .+/?o 2 e.)]. (86) 

Note that, for a propeller operating at an angle of attack, v n and / depend 
on r. rij and e* also depend on r even if the propeller is operating at zero 
angle of attack. Note also that, the definition of the source amplitude here 
differs from the standard one in that it includes the fall-off with distance 
R. 

Once source amplitude distributions Qmonopoie and Qdipoie are estimated, 
it remains to carry out the integrations in Eq. (7). Ordinarily, in the con- 
text of a frequency-domain analysis, the integral over r is computed first. 
Since, for general geometries and source amplitudes, this integral is not 
tractable analytically, it is usually computed asymptotically for, say, neax- 
or far-field observer locations for which the integrand may be simplified sig- 
nificantly. The result is given in terms of the appropriate Bessel functions. 
The remaining surface integral is then carried out analytically or numeri- 
cally depending on how faithfully the blade geometry is to be represented. 

There is, however, a less restrictive, and at the same time, more accurate 
way of computing p m B(x). This is accomplished by evaluating the r integral 
asymptotically for large-blade-count, B (or more precisely mB), using a 
modification to the standard steepest descent and saddle point methods. 
The advantage of this approach lies in the fact that it does not require 
a near- or far-field approximation and, therefore, is applicable to general 
geometries and source amplitudes. The details of the derivation will be 
presented for the zero angle of attack case first. The results will then be used 
to derive an extension to the angle of attack cases most often encountered 
in practice. For the sake of brevity, the formula will be given in terms of a 
generic source amplitude Q. The resulting expression is therefore applicable 
to both the monopole and dipole sources. 


Qmonopoie — 


Bnp 0 v n (T) 

8;r 2 


imBQ 

kR 


Qdipoie 


BQf(r) 
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1.2.1 Zero Angle of Attack Case (a = 0) 

Since, it is desirable to represent the blade geometry as accurately as pos- 
sible, the surface integrals cannot be carried out analytically and must, 
instead, be performed using a suitable quadrature scheme. For practical 
considerations, it is more convenient to evaluate the surface integral first. 
To do so, we begin by dividing the blade surface into N p small panels. If 
the typical panel size is sufficiently small the integrand may be assumed to 
be constant over the panel and, thus, the surface integral may be approxi- 
mated by some appropriate “mean value” of the integrand times the panel 
surface area AS. If the mean value is chosen to be the value of the inte- 
grand at the geometric center of the panel, designated ys, the error in this 
approximation is 0(L 2 ) for a given N p , where L = |y — y s | ma x- Of course, 
the error could be made arbitrarily small by choosing a large enough N p . 
Therefore, P m s(x) may be written as: 


p mB (x)~ Y, A S„, 

(9a) 

n p = l 


,2* 


In , = / Q(0) e mB *W do, 
J 0 

(96) 

$(0) = i ^0 + aV 1 — bcosOj , 

(9c) 

M tip _ 2 rr, 

A), \/x 2 + r 7 + r? ’ X 2 +r 2 + r, 3 ’ 

(9d) 


* c = -i, 

P0\ 

x 2 ~ 4~(*i - y*.) 2 > i 

p o> 


where for brevity and (y^ , y» 2) y a3 ) have been replaced by their 

cylindrical polar counterparts (xi,r, <j>) and (ya^r,,^,), respectively. Fur- 
thermore, the integration variable r is replaced by a new variable 6 = 
fir + <j) s — <j>. M t i p is the tip rotational Mach number (i.e., R t i p Q/Co) J 
A S nf the surface area of the n p th panel and is a convective phase fac- 
tor representing the collection of phase terms which do not depend on r. 
Q(0) is the source amplitude function written in terms of the variable 0. In 
order to simplify the notation, the explicit dependence of various variables 
on y s , the centroid of the n p th panel, is suppressed. <£(0) is the canonical 
phase function for a propeller operating at zero angle of attack, with a 
and b representing a combination of geometric, con vective and kin ematic 
factors. In deriving the above expressions, both the observer and source 
spatial coordinates are non-dimension alized with respect to the propeller 
tip radius Ru p . (Note that, Mqi = M 0 for zero angle of attack case.) 
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For a typical propfan B — 8, thus, even at the blade-passing-frequency 
(BPF), i.e., m— 1, the integrand in Eq. (9b) is highly oscillatory. Of course, 
one way to tackle the problem is to integrate the ^-integral using a suit- 
able quadrature scheme. That was indeed the approach taken in Nallasamy 
et al. [12]. Unfortunately, accurate computation of higher harmonics (i.e., 
m > 2) requires an ever increasing quadrature resolution to capture the os- 
cillatory nature of the integrand, and that may prove to be computationally 
expensive. It turns out, however, that it is possible to evaluate the integral 
I np rather accurately using the steepest descent and saddle point methods. 
Naturally, in that case, the variable 9 must be allowed to be complex and 
the integration path (i.e., [0, 2tt]) be deformed into an appropriate contour. 

Since in most practical applications of interest, propfans operate at a 
supersonic relative tip Mach number (i.e., M^ rc ^u p = (M? ip + Af^) 1 ^ 2 ), in 
the subsequent analysis it is assumed that A/( re f)ttp > 1* Replacing 9 with 
the variable, v = 9 + icr, leads to a phase function, $(*/), which is now 
also complex. It is fairly straightforward to find the saddle points of $(i/) 
(i.e., $ (u) = 0, where (') denotes differentiation with respect to v) and the 
appropriate steepest descent contours. The easiest way to find the saddle 
points is to define an auxiliary variable £ = cos i/, and rewrite $ (i/) = 0 in 
terms of 


(a&/2) 2 e - 6£ + 1 - (a6/2) 2 = 0, (10) 

where a and 6 were defined in Eq. (9d). Eq. (10) is clearly quadratic in £ 
and may readily be solved. Thus, for given observer and source locations, 
$(i/) has, in general, two simple saddle points in the interval [0,27r] which 
depending on whether the component of the source relative Mach number 
in the direction of the observer, i.e., M( r<s /).ej, is subsonic or supersonic, 
have different forms. The two are a complex conjugate pair if A/( rc ^.ej < 1, 
and are real if A/(ref)i e < > When Af( r< /).e,* = 1 (i.e., the “sonic condi- 
tion”), these two saddle points merge to give rise to a single second order 
saddle point. For a subsonic source, it turns out, that only one of the two 
saddle points lies on the appropriate steepest descent path and, hence, only 
that saddle point contributes to the integral. It should be noted that, the 
integrand also has an infinite number of branch points, coinciding with 
those of $(i/), of which only four lie in the region of interest. A judicious 
choice of the branch cuts guarantees that the contributions to the contour 
integral along the cuts in the i/-plane is exponentially small compared with 
the contributions from the neighborhood of the saddle points. The location 
of the appropriate saddle points, branch cuts and steepest descent contours 
for typical source locations are shown in Fig. (2). Note that, the choice of 
the steepest descent path depends on the type of the source under con- 
sideration. Also shown in this Figure, are the auxiliary (descent) contours 
needed for deforming the original contour, [0 , 2 7r] , into the steepest descent 
contour(s). Due to the periodicity of the integrand in 0, the contributions 
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from these contours cancel each other out exactly. 

Naturally, the asymptotic structure of I np for a sonic source point is 
different from that for a non-sonic one. In fact, it is fairly easy to show 
that the asymptotic expansion is given in terms of the inverse fractional 
1/3-powers of parameter mB for sonic sources, and in terms of the inverse 
fractional 1/2-powers of mB for non-sonic sources. Since, there is no simple 
way of constructing a composite expansion from these two expansions, to 
allow for a smooth transition through the sonic condition, they are not 
very convenient to use. Fortunately, employing a theorem due to Chester, 
Friedman and Ursell [13], it is possible to develop a uniform expansion. 
The details of the method may be found in Bleistein and Handelsman [14] 
among others. The basic idea is, however, rather simple. The phase function 
$(i/) is mapped conformally {v — ► C) into a much simpler function (i.e., a 
cubic polynomial in this case) which displays the relevant features of the 
original phase function. The region of interest in the complex i/-plane (i.e., 
the region containing the saddle points and steepest descent contours) is 
correspondingly mapped into a region in the complex £-plane. With the 
introduction of a few necessary parameters, the standard steepest descent 
and saddle point methods could then be applied to develop the required 
uniform expansion. The key definitions and parameters along with the final 
result are summarized below. 

The cubic is given by: 

= - ( J - 7 2 c) + /*, (Ha) 

#»=* J [*(* + ) + ♦(*'")]. t 5 =|N i,+ )- $ ( , '1i (ii&) 

where v + and denote the locations of the saddle points of $ in the 
complex i^-plane and [i and 7 are parameters defining the conformal map. 
Note that, 7, as given above, can take on three possible values or branches. 
The aforementioned theorem guarantees that one of the branches defines 
the desired conformal map (see Ref. [14] for details). That branch turns 
out to be the one for which y 2 is purely real for the problem at hand. 
With these parameters determined, the map may be constructed and the 
uniform asymptotic expansion carried out. In principle, the expansion of the 
integral 7„ p could be carried out to an arbitrary order in the parameter mfl, 
however, it turns out that, for most applications, the first term provides a 
very reasonable approximation for even the BPF component (i.e., m = 1). 
Therefore, in the subsequent analysis only the first term is considered. 

After a fair amount of algebra, the final result can be written as the 
following formula: 


I„ r ~ 



A[(mB) 2 / 3 7 2 ] 

(mB) 1 / 3 


+ d\ 


4[(mi?) 2 /y] V 
(mB ) 2 / 3 J ’ 


(12a) 
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( 126 ) 


j _ r 0 (7) + r 0 (-7) 

«o = o 


r 0 (7)-r 0 (-7) 

d l - nZ. 


r 0 (C) = Q(KO) 


du _ 7 2 — £ 2 

dc = ¥ 7 M0)’ 


( 12 c) 


where A,- and A] are the Airy function and its derivative, respectively, 
and do and di are coefficients in the asymptotic expansion. It is worth 
mentioning that, ±7 turn out to be the locations of the saddle points 
in the C-plane. The Airy function and its derivative provide a smooth 
transition from the subsonic portion of the blade to the supersonic portion. 
For 7 “ 0 (i.e., sonic sources) A, and A^ are 0(1) and, consequently, I Up 
is proportional to the inverse fractional 1/3-powers of mS. For 7^0 (i.e., 
subsonic/supersonic sources), and large mB ) A* is 0((mB)” 1 / 6 ) and A/ 
is 0((mB)^ 6 ) and, consequently, I Up is proportional to the inverse 1/2- 
power of mB as expected. Upon inserting for I np in Eq. (9a) from Eq. 
(12a) and adding the contributions from all the N p panels, the complex 
Fourier harmonic component, p m B( x )> can be calculated. 

At this point it should be pointed out that the analysis just described dif- 
fers from that presented by Parry and Crighton [15] in that they use asymp- 
totic forms of the Bessel functions in conjunction with the two-variable 
stationary phase method to asymptotically evaluate the surface integrals 
in Eq. (7). The key difference between the two analyses lies in their dif- 
fering points of view. While, the approach described in Ref. [15] seeks to 
take advantage of the interference between the acoustic pressure distur- 
bances arriving at the observer location emitted from the various sources 
distributed over the blade surface, the present method is built upon the 
idea that there is “self-interference” between the acoustic pressure signals 
emitted by the same source over one period of its revolution. From a prac- 
tical stand point, the goal in Ref. [15] is to circumvent the time-consuming 
computations required in a full scale analysis by developing simple approx- 
imations to the noise field of a many-bladed propeller in order to gain 
insight into the importance of various physical parameters which influence 
the noise radiation characteristics of a propeller. The aim of the current 
approach is to provide an efficient way of carrying out a full scale analysis 
which, at the same time, is also amenable to systematic parametric studies 
of the sort envisioned in Ref. [15]. 

In order to assess the accuracy of the formula given by Eq. (12a), its 
predictions for a test problem were compared with those obtained through 
direct numerical integration of the 0 integral in Eq. (9b). The test case 
chosen was an SR7 propfan operating in cruise conditions. An estimate 
for the aerodynamic loading distribution on the propfan blades, found by 
Nallasamy and Groeneweg [16] using an Euler CFD code developed by 
Whitfield et al. [17], was used along with the detailed geometric description 
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of the blade planform as input to the acoustic calculations. The results are 
summarized in Figs. (3) and (4). 

In Fig. (3), sideline directivity of the BPF sound pressure level (SPL) 
is shown for two representative sets of observer locations; a near-field set 
and a far-field set. These correspond to two different sideline distances; one 
at x 2 = 1.5 R t ip and the other at x 2 — 5 Ru p . In each case, the observer 
axial location was varied between x\ — —Rtip (i.e. , forward of the plane 
of rotation) and x\ — -j-2 R t i p (i.e., aft of the plane of rotation). The solid 
lines represent the asymptotic estimates of the SPLs, and the dotted lines 
their numerically computed values. As it can be seen, both in the near- and 
far- fields, the agreement between the numerical results and the asymptotic 
estimates is quite good, with the maximum error of about 1.5 dB occurring 
at 0.5/?tt P aft of the plane of rotation for the near-field observer. 

Due to the nature of the asymptotic expansion, the maximum error in 
approximation occurs for the lowest value of the mB parameter, which 
is m = 1 corresponding to the BPF. As m increases, the error in the 
asymptotic approximations diminishes rapidly. Therefore, it is no surprise 
that the agreement between the estimates obtained from the asymptotic 
formula and the numerically computed results improves with increasing 
harmonic order. This result is highlighted in Fig. (4), where the waveforms 
of the radiated sound at the 0.5i?*,p aft location for the two sideline dis- 
tances mentioned earlier are presented. Note that, at the axial location 
chosen, the asymptotic estimate of the BPF sound pressure level deviated 
the most from their numerically computed values (see Fig. (3)). Once again, 
the agreement is exceptionally good, particularly at the far-field location 
where the asymptotic and numerical results are virtually indistinguishable 
from each other. In constructing the waveforms the first ten harmonics were 
used, which meant that, in carrying out the numerical integration, a pro- 
gressively finer integration step size had to be utilized. This of course was 
reflected in an increase in the required computation time. The asymptotic 
estimates on the other hand, required virtually no additional time as com- 
pared with that required to compute the BPF results. This is simply due to 
the fact that the location of the saddle points and steepest descent paths, 
and the structure of the conformal map are independent of the parameter 
mB. It is worth mentioning that the amount of time required to construct 
the waveforms using the asymptotic method was more than_ an order of 
magnitude smaller than that required in the full numerical integration! 

Having demonstrated the utility of present approach, we next outline the 
extension to include the angle of attack effect. Once again, for brevity, only 
important steps in the analysis are presented. 

1 . 2.2 Non-Zero Angle of Attack Case (a^O) 

In general, the development of the analysis in this case follows closely that 
presented in the previous section. However, the details of the analysis itself 
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are somewhat more complicated. For example, the phase function $(0,a), 
(which is now more complicated than that for $(0 , 0) given by Eq. (9c)), has 
four saddle points (i.e., $ ( a) = 0 is a quartic equation). While, in princi- 
ple, it is possible to find the roots of a quartic analytically, the expressions 
involved are rather lengthy. Fortunately, since in most practical applica- 
tions of interest, the propeller angle of attack is small (typically no more 
than 5°), it is possible to find analytic approximations to the saddle points. 
This is done by treating the angle of attack case as a small perturbation 
to the zero angle of attack case (henceforth called the unperturbed case). 
The details of this development are given below where for convenience the 
analysis is carried out in terms of the parameter A/ 0 3 = Mq sin a, instead 
of or itself. 

We begin by noting that for a small angle of attack, A/03 is a small 
parameter, i.e., | Afo3 1 "Cl)* Using perturbation analysis (see, for example, 
Bender and Orszag [18]) it can be shown that two of the four roots of 
$ (v } M03) — 0, denoted by lie near those of the unperturbed problem 
^ (1/, 0) = 0 (i.e., v ± ) } while the other two scale with the inverse powers of 
M03 and, thus, lie much farther away. Therefore, since the perturbations in 
the integration contours are also small, the asymptotic behavior of I Up is 
controlled by the pair of saddle points . Of course, they could also merge 
to give rise to a second order saddle point as in the unperturbed case. 

To find i >± , first the phase function $(t/, A/03) is expanded in a Taylor 
series about A/03 — 0 as follows; 


<*(", M 03 ) = *o(*0 + M 0 3 *i(") + + 0(Mq 3 ), (13a) 


where 


0 ), 


<J 2 (^) = 


$1(1^) 


M03) 


dMo 3 


J — 0 


0Mq 3 . Af O 3 = 0 


(136) 


Qq(v) is given by Eq. (9c) (with 8 replaced by u) ) and ^x(i^) and $2 (y) can 
be easily found by carrying out the differentiations with respect to A/03* 
Since, the salient features of the analysis which follows could be demon- 
strated without the need for the explicit expressions for $1(1/) and 
they are not given here in order to conserve space. The reason for retaining 
second order terms in Eq. (13a) will become clear shortly. 

Next, Eq. (13a) is differentiated with respect to v and set equal to zero; 


A/ 03 ) ~ *;,(«/) + Afos^iW + 2 M 03^( v ) = °> 


(14) 


II 


where the higher order terms in A /03 are neglected to maintain consistency. 
Rewriting the resulting equation in terms of the auxiliary variable £ = cos v 
and employing the notation used in Eq. (10), we find that, 

(ab/2) 7 ^ - b£ + 1 - (ab/ 2) 2 + 2 M 03 [(1 - bt)Q '^coe" 1 £)] + 

M»(l - &£) [*” ( cos_1 0 + ^(c 08 " 1 0] = 0, (15) 

where, as before, only the second order terms are retained. Note that, 
when A /03 = 0, Eq. (15) reduces to the unperturbed case given by Eq. 
(10). Therefore, as pointed out earlier, it is reasonable to expect that two 
of its roots lie near those of the unperturbed case. In order to solve Eq. (15) 
for the desired roots, = cos the functions inside the curly brackets 
are expanded in a Taylor series about the unperturbed roots £* = cos . 
After collecting terms we find that Eq. (15) can be approximated by: 

[(afc/ 2) 2 + M 03 H?]S 2 - [fc - 2 M m (H? + f±J^)]f+ 

[l - (aft/ 2) 2 + Mo 3 (2Hq - 2t*H? + ^ Hf) + = 0 , (16a) 

where 

^o ± = (l-^ ± )$i(cos- 1 4 ± ), (166) 




' 1 - 6^ 

(l_^)l/2 


^(cos 1 £ ± ) + 6$ 1 (cos 1 £ 




(16c) 




^7 (cos 1 ^ ± ) 




~ ^(l-6^) 

(l_^±>)3/2 


26 

(l_£±’)l/2 



(16 d) 


K o = (! - b ^ ± ) [ $ r ( cos 1 (*) + ^ 2 ( cos_1 £*)] ■ (16c) 

Hq, Hfi Hy denote the coefficients of successive terms in the Taylor series 
expansion of the function inside the first curly bracket in Eq. (15), and 
Kq denotes the leading term in the Taylor series expansion of the function 
inside the second curly bracket. Note that, cos” 1 in Eqs. (16a-16e). 

The roots of this “new” quadratic can now be easily found. Since, Eq. 
(16a) has two roots for each of and we choose the root that reduces 
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to £ + or in the limit of — ► 0. As expected, turn out to be only 

slightly different from the unperturbed roots It is worth mentioning that 
alternative perturbation methods could be used to find However, the 
particular approach described above provides a uniform representation of 
the roots even when the unperturbed part of Eq. (16a) has a single double 
root (i.e., the sonic condition for the unperturbed problem). Near that 
condition the accuracy of the approximation drops to 0(Mq g 2 ). This is why 
in the original expansion (Eq. (13a)) the second order terms were retained. 
If only the first order terms were kept, near the sonic condition, the accuracy 
would be no more than With the two perturbed saddle points 

now known, the asymptotic expansion of can be carried out following 
the same procedure as was outlined in the previous section; the difference 
being that $(*/, a) is now replaced by the three terms on the right hand 
side of the Eq. (13a), and are replaced by . 

Preliminary applications of the angle of attack formula to a few test cases 
have been rather encouraging. As in the case of zero angle of attack, it pro- 
vides accurate estimates for the integral 7„ p . The utility of the method in 
solving realistic problems is demonstrated with the following example which 
involves the prediction of the spectrum of noise of SR7 propfan operating 
in cruise and at an angle of attack of a = 1.6°. The unsteady aerodynamic 
loading distribution for this case was computed numerically by Nallasamy 
and Groeneweg [19] and was shown (see Heidelberg and Nallasamy [20]) 
to be in quite reasonable qualitative , as well as quantitative, agreement 
with the experimentally measured airloads. Using this aerodynamic load- 
ing distribution along with a detailed geometric description of the blade 
geometry as input, calculations for the spectrum of the noise of the prop- 
fan were carried out utilizing the asymptotic formula. The results were 
then corrected for non-free-field effects (such as scattering) which are not 
accounted for in the theoretical model, The corrections used, which were 
found in Spence [21], were 1.0, 3.6 and 5.3 dB for the first three harmonics, 
respectively. For higher harmonics a uniform correction of 5.0 dB was as- 
sumed. The experimentally measured spectrum corresponding to this case 
was presented in Ref. [22], as part of the Propfan Test Assessment (PTA) 
program (see Poland et al. [23]). The comparison between the measured 
and asymptotically computed spectra is presented in Fig. (5). In general, 
the agreement is quite good with appreciable deviations occurring only at 
higher harmonics where the neglected quadrupole sources probably play an 
important role. 


1.3 Concluding Remarks 

An asymptotic approach was presented which allows for accurate and ef- 
ficient calculation of the noise field of a propfan at angle of attack with- 
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out the need for simplifying assumptions usually employed in propeller 
noise theories. A closed form expression, involving the Airy function and 
its derivative, gives a uniform representation of the propeller noise field 
characteristics in both the near- and far-fields. The preliminary results 
shown here indicate that the method provides a useful theoretical tool in 
propeller noise research. 
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Fig. (1) Definitions of Geometry and Coordinate System. 
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Fig. (2) Saddle Points and Steepest Descent Contours in the Complex 
j/-Plane. 
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Fig. (3) Comparisons of Predicted Sideline Directivities of SR7 Propfan 
at Cruise Conditions and Zero Angle of Attack (a = 0). 
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Time, milliseconds 


Fig. (4) Comparisons of Predicted Pressure Waveforms for SR7 Propfan 
at Cruise Conditions and Zero Angle of Attack (a = 0). 
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